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We summarize recent results concerning the evolution of second order perturbations 
in flat dust irrotational FLRW models with A / 0. We show that asymptotically these 
perturbations tend to constants in time, in agreement with the cosmic no-hair conjecture. 
We solve numerically the second order scalar perturbation equation, and very briefly 
discuss its all time behaviour and some possible implications for the structure formation. 
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1. Introduction 

First order perturbations of Friedmann-Lemaitre- Robertson- Walker (FLRW) mod- 
els fail to account for a number important relativistic properties that arise when 
non-linearities are taken into account. An important example of this is the so called 
mode coupling. Even if one considers purely scalar perturbations, a comoving scale 
k at second order is sourced by any other scale k' . In addition, the scalar, vector 
and tensor perturbations which are decoupled to the first order, become coupled 
once nonlinear perturbations are taken into account. Thus for example taking into 
account such couplings at second order results in the initial pure scalar metric 
perturbations to necessarily generate second order vector and tensor perturbation 
modes (see e.g.El). In view of this, it is important to develop higher order perturba- 
tion schemes which go beyond the first order and can thus account for the physical 
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effects which are not taken into account up to the linear order. This would also 
improve the level of accuracy of results as well allowing the study of stability of the 
results obtained using first order perturbations. 

Among concrete motivations for the use of a second order perturbative scheme 
is the need for more accurate results from the new generation of gravitational wave 
detectors as well as the increased computational precision required in order to 
analyse the cosmic microwave background anisotropics .cm Both these problems are 
likely, in future, to require analyses that go beyond the linear order. In addition, the 
non-linear analysis in the context of cosmoiogicai modeling is also relevant on scales 
much smaller than the cosmoiogicai horizon where it might substantially modify 
the first order resultsQ and therefore, as we shall also see below, become important 
for the study of the structure formation. Furthermore, as far as observations are 
concerned, second order effects can be important on small scales in order to account 
for non-linear effects such as the Rees-SciamaO and gravitational lensing effects Isee 
e.g.B). For a list of references concerning these questions we refer the reader tcO. 

Second order perturbations have so far not been widely studied in general rela- 
tivity. In the cosmoiogicai context they seem to have been first employed by Tomitai 
to study the evolution of scalar perturbations in the Einstein-de-Sitter model using 
a synchronous gauge. This study was repeated by Mataresse et. al.i, who obtained 
similar results using a comoving approach and Russ et al.i who have included sec- 
ond order terms resulting from a coupling between groszing and decaying scalar 
perturbation modes. There have also been recent worka!3 showing that the back 
reaction of the second order perturbations (to be precise, perturbations quadratic 
in first order terms) are important in early universe scenarios. Furthermore, it has 
been shown that in inflationary models the magnitude of the second order pertur- 
bations can be comparable to those Jm the first order and therefore a non- linear 
perturbative analysis may be crucialO 

The domain of applicability of the theory of second order perturbations that we 
shall describe is that of small perturbations about an homogeneous and isotropic 
background. We shall employ it as an approximate framework in order to study 
the effect of non-linearities in structure formation scenarios as well as the study of 
the stability of cosmoiogicai results obtained using first order perturbative schemes 
with respect to nonlinear perturbations. 

In this paper we shall summarise recent results concerning the exact asymptotic 
solutions of the second order perturbation equations. An immediate consequence 
of these results was to prove the non-linear asymptotic stability of the de-Sitter 
attractorE3, and thus generalise the cosmic no-hair conjecture. In general, how- 
ever, the perturbation equations cannot be solved analytically. We shall show pre- 
liminary results concerning the numerical integration of the scalar perturbation 
equation for all times. Finally, we discuss possible implications on the structure 
formation. We use units in which 8irG = c = 1. Greek indices take values 1, 2, 3. 
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2. Perturbation Equations 

We consider a homogeneous and isotropic background space-time given by a flat 
FLRW metric with an irrotational dust source field. The relevant equations for the 
expansion 9 a p in presence of the cosmological constant A are given by the evolution 
equation, the Raychaudhuri equation and the energy constraint: 

o a p + ee aP + R* a p = -p5 aP + As a p (i) 
+ o aP e a p + \p = a (2) 

2 -9 af3 6 af3 +R* = 2(p + A) (3) 

where the dot denotes a derivative with respect to coordinate time t, p is the matter 
density and R* is the 3-Ricci scalar. We use the notation of Matarrese, Mollerach 
and BruniEJ and generalise their formalism to the case with nonzero A. 
Using a conformal rescaling we can write the metric as 

ds 2 = a 2 (r)(-dr 2 + j a/3 dx a dx ), 

where a is the scale factor. We work in the synchronous gauge (which for dust is 
also comoving) in which the first order metric perturbations can be written as 

7a/3 = <W + 7«J with = ~2(j)5 a p + D a pX + (4) 

where D a p — d a dp — ^S a pV 2 , the superscript W denotes first order quantities, <fi 
and x are the trace and tracefree parts of the scalar perturbation modes and the 
tensor modes are represented by TT a p 7 which is transverse and divergence free. In 
the linear theory, the scalar, vector and tensor modes are independent and for an 
irrotational space-time one can set the vector modes equal to zero. 

We recall that the evolution equation for the tensor modes can be obtained from 
the tracefree part of (|l]) in the form 

<P + ^' a p ~ V 2 TT af3 = 0, (5) 

where the prime denotes a derivative with respect to conformal time r. Using tensor 
harmonics, which are eigenfunctions of the Laplace operator, one can set V 2 TT a p = 
—k 2 Tr a p, where k is the order of the harmonic. 

The first order scalar perturbation equation can be obtained from the continuity 
equation (which in our case reads p + Bp = 0), together with ([!]) and (g), to give 

5" + - \a 2 Pb 5 = 0, (6) 
a 2 

where 5 = (p — Pb)/pt and p b denotes the background matter density. One can 
further make use of the freedom of the synchronous coordinate system to specify 
S = — iV 2 %. This then completely fixes the gauge. 
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Having recalled the first order perturbation equations we shall now consider 
the second order case. The conformal spatial metric tensor of the second order 
perturbed space-time can be written, as 

7a/J = *a/s+7$ + 57S! 72 = -2^ (2) <5 Q/3 + xi% (7) 

where (f>^ is the second order scalar perturbations and x„p ^ s trace-free. In order 
to obtain a simpler form for the second order perturbation equations we shall not, 

(2) 

as in the first order case, separate the trace-free part of the \ a p perturbations. 
Substituting (jjj) in (0) we obtainEB the following evolution equation for <j)^ 



+ Pba 



If (Da (l)<x\ , If (Dad (!) (Dap (1) \ r / [1) (lk 
J I 7^° " 7oa J + 2 ( 7 l j ^7a/3 - 7o 7oa/3 J - ^0 [li > - 7oa 
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where the subscript "0" denotes the evaluation at an initial time t = tq. 
3. Exact Asymptotic Results 

It is well known that the first orderoerturbations of flat dust FLRW models asymp- 
totically tend to constants in time.Ej Now, in order to test the stability of this result 
with respect to second order perturbations we can solve the asymptotic form of (||) . 
Keeping the lowest order terms on the right hand side of Eq. (||) we have 

W" _ V)' + l^por^ = F(x)r + 0(r 2 ), (9) 



.3/2 



(10) 



where F is a spatial function. The solutions to (||) have the form 
0( 2 )(r,x) = ^M +r C 3 (x)J (f,^T 3 / 2 ) +rC 4 (x)F Q,|at 

where C3 and C4 are arbitrary functions and J and Y are Bessel functions. There- 
fore, asymptotically^ approaches a constant value in time. 

The evolution equation for x^p can be obtained from (|l|) (seetil) and its asymp- 
totic form (for each a and 0) is given by 

X$" - 2 ~xS + fc\S = E ^k) + 0(r), (11) 
where i£ Q| a depends on fc. This equation can be solved to give 

X a 2 p( T 'k) = Eafi }^ > + c a p(kTCOs(kT) -sin(fcr)) + d a p(kT sin(fcr) +cos(fcr)), (12) 



a Note that t — > 00 corresponds to r — > 0. 
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where c a p and d a fj are arbitrary. As a result, as r — > 0, Xq,^ tends to a constant in 
time which depends on the asymptotic values of (p^ and c£ Q ^. On the other hand, 
these values are related to the first order initial free data so, at second order there 
are four independent functions. 

To summarise, we showed that second order scalar and tensor perturbations 
asymptotically approach a constant value in time. This demonstrates that the 
asymptotic behaviour of the first order perturbations is stable to the presence of sec- 
ond order perturbations. In addition, using an invariant analysis (seec3; cf. alsoEj), 
we established the validity of the cosmic no-hair conjecture in nonlinear settings. 

4. Numerical Results 

Our analysis so far has concentrated on the asymptotic behaviour of the pertur- 
bations. The all time evolution of the second order perturbations is of potential 
importance for processes that evolve in time, such as the structure formation pro- 
cess. Here as a first step in this direction we numerically integrated the equation 
(||) for several initial conditions in order to investigate the evolution of in time. 
We found that even though asymptotically there is a single future attractor that 
determines the qualitative behaviour of (f)^ at late times, the precise approach 
to this attractor can be different at intermediate times, as shown in Figure [l]. 
The two curves represent the evolution of for different initial conditions. Both 
curves have initial values A = 0.01, a(r ) = 0.9, (2) (to) = 0.1, (1) (to) = 0.01, 
4^ (to) = 0.01 and po = 0.01. The lower curve has <^ 2 ' (tq) = 0.01, and the upper 
curve 0'( 2 )(t o ) = 0.025. Varying <^' (2) (to) between 0.1 and 0.025 results in a family 
of curves lying between the two depicted curves. As can be seen, the evolution of 
0( 2 ) can have distinct behaviours over intermediate times; having for example either 
one or three inflection points depending on the initial data. 

This is an example of a dynamical system with the same attractor but different 
transient phases. This type of behaviour could result in different rates of structure 
formation over intermediate time scales, depending upon the choice of the initial 
data. In particular we could have a possible intermediate period where <f>^ can be 
almost constant in time, as is the case with the lower curve in Figure [I]. 

5. Discussion and Conclusions 

We have summarised our very recent resultsElQ on the evolution of second order 
perturbations in flat irrotational dust FLRW models. We have shown that second 
order perturbations tend asymptotically to constants in agreement with the cosmic 
censorship conjecture. 

In order to obtain the behaviour of such perturbations at intermediate times, 
we have also integrated numerically the scalar perturbation equation and found 
different transient behaviours over intermediate time scales, depending upon the 
choice of the initial conditions. This can have interesting consequences for nonlinear 
phenomena that evolve in time, such as structure formation. 
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Fig. 1. Phase portrait representing the evolution of second order scalar perturbations <j>^ in a 
comoving gauge for a flat FLRW dust model with A = 0.01. 
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